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Abstract 

Given  a  time  history  of  desired  moments,  the  control  allocation  prob¬ 
lem  is  to  solve  for  the  effector  inputs  so  that  some  norm  of  the  error  between 
the  achieved  and  desired  moments  is  minimized.  Existing  methods  solve  for 
the  actuator  deflections,  while  accounting  for  Magnitude  and  rate  limitations 
of  the  effectors.  In  this  paper,  we  propose  the  Dynamic  Control  Allocation 
(DCA)  Method,  that  also  accounts  for  effector  dynamics,  in  addition  to  mag¬ 
nitude  and  rate  limits.  We  show  through  numerical  experiments  that  the  DCA 
method  allocates  the  desired  moments  according  to  effector  bandwidths  -  that 
is  the  slow  effectors  are  allocated  the  lower  frequencies  in  the  desired  mo¬ 
ments.  The  numerical  simulations  also  show  that  the  DCA  outperforms  the 
existing  simplex  algorithm  based  LP  method,  that  does  not  account  for  actua¬ 
tor  dynamics. 
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Nomenclature 

Dynamic  Control  Allocation 
3-vector  of  desired  moments  at  time  t 
3-vector  of  moments  achieved  by 
the  effectors  at  time  t 
Moment  error  =  Gs,des{t)  ~ 

The  p  norm  of  the  function  /  for  1  <  p  <  oo 
The  Fourier  transform  of  the  function  / 

The  space  of  functions  defined  on  the 
interval  I  with  bounded  p  norm 
The  functional  to  be  optimized  as  a 
function  of  the  control  input  u 

1.  Introduction 

The  problems  of  control  allocation  and  reconfigurable 
control  have  recently  been  widely  studied.  A  review 
of  existing  methods  can  be  found  in  Bodson  [1].  Due 
to  onboard,  real-time  computational  constraints,  exist¬ 
ing  methods  address  the  problem  of  control  allocation  at 
each  discretized  time-step  separately  and  attempt  to  min¬ 
imize  the  difference  between  the  desired  and  the  achiev¬ 
able  moments,  while  accounting  for  rate  and  magnitude 
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limits  on  the  effectors.  This  approach  can  be  found  in 
Buffington  [2],  and  Doman,  Ngo  [3].  The  downside  of 
this  approach  is  that  the  dynamics  of  the  actuators  are 
not  taken  into  account,  and  this  could  mean  that  the  mo¬ 
ments  actually  achieved  might  be  significantly  different 
from  the  computed  moments.  This  becomes  a  problem 
in  reconfigurable  control  because  the  goal  of  reconfigura¬ 
tion  is  to  recover  from  damaged  effectors.  Burken  et  al. 
[4]  and  Pachter  et  al.  [5]  try  to  include  actuator  dynam¬ 
ics  by  first  solving  an  LQR  problem  and  then  solving  for 
the  closest  achievable  moments  when  rate  and  magnitude 
limits  are  present.  We  propose  to  solve  for  the  inputs  to 
the  effectors  in  one  step  via  an  optimization  procedure. 

In  this  paper,  we  propose  the  Dynamic  Control  Allocation 
(DCA)  method  that  takes  into  account  individual  effector 
dynamics  as  well  as  rate  and  magnitude  limits.  We  uti¬ 
lize  an  effector  model  that  incorporates  these  effects,  to 
predict  the  moments  achieved  for  some  input.  The  DCA 
method  minimizes  the  difference  between  the  desired  and 
’’predicted”  achievable  moments  over  all  possible  inputs 
to  the  effector  model.  For  effectors  whose  dynamics  can 
be  modeled  by  a  linear  system,  our  method  leads  to  a 
convex  optimization  problem. 

A  different  approach  was  taken  in  our  earlier  work  [6] 
where  we  assumed  the  control  allocation  was  solved  by 
existing  methods  and  constructed  a  controller  for  an  effec¬ 
tor  with  magnitude-limited  first-order  dynamics  so  that 
it  can  follow  the  commanded  deflection.  However,  the 
analysis  gets  quite  complicated  for  effectors  with  higher- 
order  dynamics,  while  the  method  presented  here  is  less 
complex. 

Mathematical  Preliminaries 

Let  I  =  [toj^i]  he  an  interval  of  time  over  which  the 
control  allocation  problem  is  required  to  be  solved.  Let 
Gs,des{t)  denote  the  3-dimensional  vector  of  desired  mo¬ 
ments  and  5(t)  denote  the  n-dimensional  vector  of  control 
effector  positions  at  time  t  €  /.  We  assume  n  >  3.  The 
space  of  functions  Gs{t)  and  Gs^desil^)  where  t  G  I,  is  a 
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Gs^desil') 

G6{t) 

E{t) 

ll/llp 
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J{u) 
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vector  space  on  which  we  can  consider  several  types  of 
norms.  The  useful  norms  that  could  be  considered  are 
defined  as  follows: 

ii/iIp  =  ymrdty  (1) 

ll/lloo  =  ess  sup \f{t)\  (2) 

t£l 

If  ll/llp  <  oo  for  some  p  such  that  1  <  p  <  oo  then  /  is 
said  to  be  in  the  normed  vector  space 

The  control  allocation  problem  can  be  formulated  in  the 
time-domain  in  terms  of  the  moment-error,  or  in  the 
frequency  domain  by  considering  the  Fourier  transform 
of  moment-error.  The  Fourier  transform  of  a  function 
/  e  L^{I)  is  the  function  /  defined  by: 

=  (3) 


2.  Control  Allocation 

Suppose  there  are  n  control  effectors.  The  effectors  are 
assumed  to  satisfy  the  differential  equations: 

=  fiMt),Ui{t)),  (11) 

Si{t)  =  hi{xi{t)yUi{t))  i  =  tel  (12) 

The  vector  functions  fi  and  hi  could  be  nonlinear  due  to 
rate  and/or  magnitude  limitations.  The  control  effective¬ 
ness  function  J5(p(t),*)  maps  5{t)  =  [Si{t)  •  ••  ,  <5n(t)]^  to 
the  3-dimensional  vector  of  moments  Gs{t)  produced  by 
the  effectors: 

Gsit)  =  B{p{t),Sm  03) 

where  p(t)  is  a  set  of  parameters  such  as  Mach  number, 
angle  of  attack  and  side  slip  angle.  It  is  assumed  in  this 
work  that  the  time  scale  on  which  B(p(*),*)  changes  is 
much  larger  than  the  length  of  the  interval  I.  This  as¬ 
sumptions  is  justified  because  the  parameters  p{t)  are 
typically  slowly  changing.  Thus  we  regard  5(p(t),*)  to 
be  denoted  by  jB(-),  so  that: 


where  j  =  y/^.  Then  we  have  the  following  useful  theo¬ 
rem: 

Theorem  LI:  (Rudin  [7]) 

•  If  /  C  L^{I),  then  /  G  Co(R),  and  ||/lloo  <  ll/lli- 

•  If  /  G  1/2(7),  then  1|/||2  =  \\fh  (Parseval-Plancheral). 

Here  Cb(R)  is  the  supremum^normed  Banach  space  of  all 
complex  continuous  functions  on  R  that  vanish  at  infin¬ 
ity.  The  Fourier  transform  of  a  vector- valued  function  is 
defined  component-wise.  If  E{t)  =  [Ei{i)  7?2(t)  £'3(t)]^, 
then  E{j^)  ==  If  Ei{t),  i  = 

1,2,3  G  L^{I)  or  7/^(7),  then  Theorem  1.1  can  be  ap¬ 
plied  as  follows: 


li^(-)lloo  := 

.m£^3  ll7(-)lloc 

(4) 

< 

.m^j  Pi(-)lli 

3 

(5) 

< 

i=l 

(6) 

= 

(7) 

11^(0112  - 

f  E  ll^i(')llij 

1  (8) 

/  ®  ^ 

h 

= 

(  E 

1  (9) 

= 

P(-)i|2. 

(10) 

In  the  control  allocation  problem,  E{t)  is  taken  to  be  the 
difference  between  the  desired  moments  and  the  achieved 
moments. 


Gsit)  =  B{S{t)y  tel.  (14) 

Note  that  this  general  definition  can  also  be  used  to  model 
actuator  interactions.  Define  E{-)  =  G,5,des(*)  “  Gs{^)  to 
be  difference  between  the  desired  and  achieved  moments. 
As  usual,  we  denote  E{u;)  as  the  Fourier  transform  of 
Eit). 


The  control  allocation  problem  can  be  posed  in  several 
ways: 

Problem  Statement  1:  Obtain  the  control  inputs  to  the 
effectors  Ui{t);  i  =  1,  •  •  -  n,  t  G  7,  so  that  |17;(‘)llp  is  min¬ 
imized  for  some  p  such  that  1  <  p  <  cx). 


Problem  Statement  2:  Obtain  the  control  inputs  to  the 
effectors  Ui(t);  i  =  1,  •  •  -  n,  7  G  7,  so  that  ||E(-)llq  is  min¬ 
imized  for  some  q  such  that  1  <  ^  <  oo. 

To  begin  the  development  of  the  DCA,  denote  Jp{u)  = 
||7?(-)||p  and  Jg{u)  =  ttJ5(-)|lg-  Also  denote  Up  = 
arg  min  Jp{u)  and  Up  =  arg  min  Jp{u),  Then  it  is  clear 
ueL^i)  ueL^ii) 

from  the  discussion  in  the  mathematical  preliminaries 
that: 

•  ll/lloo  <  Pill  ==^  mm  7oo(^)  <  nun  Ji(w); 
•\\M\2-\\E\\2  irnn  J2(u)  =  imn  J2(u), 

From  the  second  result,  we  see  that  the  same  solution  U2 
minimizes  the  2-norm  of  the  moment-error  or  its  Fourier 
transform.  The  same  is  true  when  min  Ji(u)  =  0  (or 

u 

equivalently,  minJoo{u)  =  0).  Constraints  on  the  set 

u 


2 


in  which  u  belongs  could  lead  to  mmJi{u)  7^  0.  By 

u 

our  notation,  ui  =  arg  min  Ji(u),  and  it  is  possible 
neLi(/) 

that  there  exists  Uoo  —  arg  min  Joo{u)  that  satisfies 

uGL^I) 

Joo{uoo)  <  onr  goal  is  to  allocate  controls  ac¬ 

cording  to  bandwidth  of  the  actuators,  then  it  seems  that 
a  more  natural  cost  function  would  be  Joo(^))  rather  than 
Ji(tt).  We  call  the  solution  Uoo  CIS  the  optimal  solution.  How¬ 
ever,  as  the  minimization  problem  is  more  easily  solved  in 
the  time-domain  (which  leads  to  the  solution  ki),  we  ob¬ 
tain  a  sub-optimal  solution.  The  1-norm  is  used  frequently 
in  Control  allocation  literature  as  it  leads  to  linear  pro¬ 
gramming  approaches  to  the  numerical  solution  [1],  [2], 

[3] , 

Existing  methods  allocate  control  at  a  time  instant  t  by 
minimizing  the  moment-error  at  that  time  alone  [1],  [2], 

[4] ,  [5],  [3].  One  can  see  that  as  the  interval  I  in  our 
discussion  collapses  to  one  point,  we  obtain  the  existing 
methods.  Some  methods  take  into  account  actuator  rate 
and  magnitude  limits  in  the  control  allocation  problem 
but  do  not  include  actuator  dynamics  [1],  [2],  [3].  Other 
methods  try  to  include  actuator  dynamics  by  first  solving 
an  LQR  problem  and  then  solving  for  the  closest  achiev¬ 
able  moments  (using  LP  or  QP)  when  rate  and  magnitude 
limits  are  considered  [4],  [5]. 


3.  Numerical  Experiments 

We  consider  control  effectors  such  as  flaperons,  rudders 
only  in  the  following  numerical  study  We  assume  the 
air  vehicle  is  on  an  unpowered  descent  so  that  the  engine 
is  not  used,  though  it  should  be  noted  that  our  method 
is  general  enough  to  handle  complex  engine  models  also. 
Each  effector  is  considered  to  be  a  2nd  order  linear  sys¬ 
tem  with  magnitude  and  rate  limitations  as  shown  in  Fig¬ 
ure  1.  In  the  figure,  the  signals  correspond  to  uj  in 
Equations  15  and  16.  The  constants  for  the  actuators  are 
shown  in  Table  1. 


The  control  effectiveness  matrix  was  taken  to  be: 


1111 
0.5  0.5  1  1 

2  2  11 


(19) 


The  numerical  experiments  were  conducted  with  to  equal 
to  the  current  time  and  the  interval  I  only  included  the 
current  time,  so  as  to  allow  a  comparison  with  existing 
methods.  The  time-axis  discretization  used  was  0.02  sec¬ 
onds.  The  minimization  of  the  cost  function  18  with  p  =  2 
was  performed  with  the  *ucSolve’  routine  in  TOMLAB. 
The  constants  K\  and  K2  in  Figure  1  are  computed  using 
the  values  for  the  natural  frequency  w  and  damping  coef¬ 
ficient  C  in  Table  1  according  to:  Ki  —  and  ^^2=0;^. 


Our  proposed  method  includes  both  actuator  dynamics 
as  well  as  rate  and  magnitude  limits  on  the  actuators. 
Suppose  that  the  Equations  11  are  integrated  to  yield: 


Xi{t) 

=  /  fi(xuUi)dt 

(15) 

Jto 

Si{t) 

•  .n  (16) 

where  t  G  /.  We  consider  the  “current”  time  to  be  in  I 
and  the  problem  is  to  solve  for  u{t)  for  the  entire  inter¬ 
val  /.  Once  this  is  done,  the  inputs  corresponding  to  the 
current  time  is  applied  to  the  actuators. 

We  discretize  the  time  axis  into  instants  r^,  fc  =  0,  •  ♦  •  , 
such  that  To  =  to  and  =ti.  The  problem  is  then: 

min  J{u)  =  min  ||£?(-)IIp  (^7) 

=  min  \\B{5{-))^GsM')\\p^  0^) 

where  p  =  1  or  2.  There  is  no  constraint  on  the  set  to 
which  u  belongs.  When  the  system  dynamics  are  linear, 
then  <5(*)  is  a  linear  function  of  u(*).  Furthermore,  if  the 
achieved  moments  B(^(‘))  can  be  approximated  by  a  lin¬ 
ear  function  Bi  (5(-)  where  is  a  3  x  n  matrix,  then  J{u) 
is  a  convex  functional  of  the  input  function  u(*).  The  lin¬ 
ear  approximation  of  the  control  effectiveness  function  is 
valid  when  /  is  a  small  enough  interval  and  can  be  seen 
in  works  of  other  researchers  [1],  [2],  [3]. 


Figure  2  shows  the  result  of  the  numerical  experiments 
with  /  =  {to}.  The  desired  moments  Gs.desi')  were  taken 
to  be  2sin(^)[l  1  1]^  -h  [3.2  2.4  4,8]^,  where  t/  =  9. 
Figure  2(a)  shows  the  desired  moments  and  the  moments 
achieved  by  our  proposed  Dynamic  Control  Allocation 
(DCA)  method.  It  also  shows  the  moments  commanded 
and  achieved  by  the  simplex  method  based  LP  algorithm 
in  Bodson  [1].  For  the  latter  method,  the  commanded 
deflections  6des  were  computed  only  considering  the  rate 
and  magnitude  limits  on  the  effectors  and  ignoring  the 
dynamics.  The  commanded  moments  in  Figure  2(a)  are 
then  given  by  BSdes-  Then  the  commanded  deflections 
are  applied  to  a  model  of  the  actuators  and  the  achieved 
moments  are  computed  according  to  B5.  One  can  see 
that  though  B  6des  tracks  the  desired  moments  well,  the 
achieved  moments  are  far  from  satisfactory.  On  the  other 
hand,  the  DCA  method  produces  achieved  moments  that 
track  the  desired  moments  better. 

Figure  2(b)  shows  the  commanded  and  response  deflec¬ 
tions  of  the  effectors.  One  can  clearly  see  that  initially  the 
slow  changing  moments  are  allocated  to  the  rudders  (with 
the  slow  dynamics).  As  time  increases,  the  rudders  prove 
incapable  of  tracking  the  faster  changing  signals  and  the 
flaps  take  up  the  slack  at  an  increasing  rate.  Thus  both 
optimal-tracking  and  allocation  of  moments  according  to  the 
bandwidth  of  the  actuators  are  achieved  simultaneously.  One 
can  also  see  that  the  commanded  signals  5des  violate  both 
the  rate  and  magnitude  limits  as  opposed  to  the  existing 
LP  method;  however,  the  actual  rate  and  position  limits 
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Figure  1 .  Effector  schematic  diagram 


Natural  Freq. 
radjsec 

Damping  Coeff. 

Rate  Limit 
Deg  !  sec 

Magnitude  Limit 
(low, high)  Deg 

Flap  1 

15 

0.75 

4 

(0.  1.5) 

Flap  2 

15 

0.75 

4 

(0,  1.5) 

Rudder  1 

5 

0.75 

2 

(0,  1.5) 

Rudder  2 

5 

0.75 

2 

*  (0.  1.5) 

Table  1.  Actuator  parameters  for  Figure  2. 


are  satisfied  by  S. 

4.  Conclusion 

In  this  paper,  we  proposed  a  new  control  allocation 
methodology  that  takes  into  account  rate  and  magnitude 
limits  on  the  effectors,  as  well  as  their  dynamics.  Exist¬ 
ing  methods  only  take  into  account  rate  and  magnitude 
limits.  We  solve  for  the  inputs  to  the  effectors  so  that  the 
achieved  moments  match  the  desired  moments  over  an  in¬ 
terval  of  time.  One  of  the  interesting  results  obtained  was 
that  of  frequency  separation  -  the  slower  effectors  were  al¬ 
located  the  lower  frequencies  in  the  desired  signal  while 
the  faster  effectors  were  allocated  the  higher  frequencies. 
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